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One-dimensional multi-component Fermi or Bose systems with strong zero-range interactions can 
be described in terms of local exchange coefficients and mapping the problem into a spin model is 
thus possible. For arbitrary external confining potentials the local exchanges are given by highly 
non-trivial geometric factors that depend solely on the geometry of the confinement through the 
single-particle eigenstates of the external potential. To obtain accurate effective Hamiltonians to 
describe such systems one needs to be able to compute these geometric factors with high precision 
which is difficult due to the computational complexity of the high-dimensional integrals involved. An 
approach using the local density approximation would therefore be a most welcome approximation 
due to its simplicity. Here we assess the accuracy of the local density approximation by going 
beyond the simple harmonic oscillator that has been the focus of previous studies and consider some 
double-wells of current experimental interest. We find that the local density approximation works 
quite well as long as the potentials resemble harmonic wells but break down for larger barriers. In 
order to explore the consequences of applying the local density approximation in a concrete setup 
we consider quantum state transfer in the effective spin models that one obtains. Here we find that 
even minute deviations in the local exchange coefficients between the exact and the local density 
approximation can induce large deviations in the fidelity of state transfer for four, five, and six 
particles. 

PACS numbers: 67.85.-d,75.10.Pq,03.67.Lx 


I. INTRODUCTION 

The field studying atomic gases at extremely low tem¬ 
peratures has seen riveting developments over the past 
decade [Ij. This has pushed the field into a rather unique 
position where cold atoms in various geometries can be 
used to do quantum simulation of widely used mod¬ 
els from other fields including condensed-matter physics 
i i, particle and high-energy physics , nuclear 

physics 0, and even chemistry Q- A remarkably use¬ 
ful feature of cold atomic gas experiments is the control 
over external trapping parameters that makes it possible 
to explore physics in different dimensionalities [l| . In par¬ 
ticular, a number of impressive experiments have studied 
various aspects of interacting bosons in one dimension 
0. Using the atomic interaction resonances caused by 
one-dimensional confinement [a it has thus become pos¬ 
sible to realize many interesting one-dimensional systems 
includin g th e famous hard-core bosonic Tonks-Girardeau 
gas [a 113 (s®® Refs. (lol . [Tll|). Most recently, one¬ 
dimensional systems of interacting fermions have been 
also realized in experiments [l8| . 

In the last few years it has become possible to con¬ 
trol the particle number in one-dimensional experiments 
with great accuracy [l^ , allowing one to build a control¬ 
lable few-body system of fermions and study fermioniza- 
tion for strong interactions [a |- pairi ng 1^1 . impurity 
physics lai) ^ two-site Hubbard model |23|. and Heisen¬ 
berg spin models [23 • These developments have sparked 
a great deal of excitement in the community studying 
few-body physics and its relation to many-body phe¬ 
nomena. Over the last decade a lot of new aspects of 


these one-dimensional systems have been covered by dif¬ 
ferent authors. This includes the physics of small trap ped 
bosonic systems (single- and two-component) (25l - [^ . de¬ 
tails of few-fermion systems with two values of an internal 
degree of freedom (spin) [H, [33, iliSl and various as- 
pects of the transition from few- to many-body physics 
[52 - 1^ . Furthermore, a number of studies have 

looked into mixed systems of bosons and fermions and 
systems with particles of unequal mass . 

In the present paper we are interested in studying 
one-dimensional fermions or bosons with strong repulsive 
short-range interactions. As has been recently discussed, 
in the strongly interacting limit the system may be de¬ 
scribed by an effective Hamiltonian which has the form 
of an (anisotropic) Heisenberg spin model [13 [13 • R 
was predicted theoretically that by tuning the interaction 
strength from the weak (repulsive) to the strongly inter¬ 
acting regime for a two-component Fermi system in one 
dimension one may arrive in the ground state of an an¬ 
tiferromagnetic Heisenberg model [13 and recently this 
was confirmed in experiments [13 • However, the Heisen¬ 
berg model obtained is non-trivial in the sense that the 
exchange couplings are determined by the local geome¬ 
try of the trapping potential. Computing the local ex¬ 
change coupling constants is a formidable numerical task 
and thus several papers have discussed the possibility to 
use various approximations to obtain these quantities. 
A very neat approximation is the strong-coupling ansatz 
of Levinsen et al. [ssl . [b^ which allows one to get an 
extremely accurate set of exchange coefficients for arbi¬ 
trary system sizes in the case where the external confine¬ 
ment is given by a harmonic oscillator. It has also been 
conjectured that for smooth potentials, the well-known 
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local density approximation (LDA) should give a very 
accurate value for the local exchange coefficients 1 ^ 
and numerical results for up to six particles show that 
the deviations are at most a few percent [s^- This is 
a very reasonable expectation and many studies in the 
past have used the LDA when studying th e prop erties of 
many-particle systems in one dimension |85l - l93j . How¬ 
ever, this does not necessarily imply that the LDA works 
equally well for smaller systems as large system sizes can 
sometimes average out some of the finer details that the 
LDA may not capture. 

In the present paper we test the performance of the 
LDA for strongly interacting one-dimensional systems 
with up to six particles by comparing it to exact calcula¬ 
tions. We do this for several different potential forms, in¬ 
cluding some experimentally relevant double-well geome¬ 
tries. To the best of our knowledge, no previous study 
has presented results for the effective spin models with up 
to six particles in non-harmonic confinement. We gauge 
the performance of the LDA against exact results not 
only for statics (producing the local exchange constants 
needed for the effective Hamiltonian) but also for dynam¬ 
ics. For the latter case we study quantum state transfer 
using the spin models that one obtains with the LDA 
and with exact calculations. Transfer of quantum states 
in two-component spin systems is a delicate process that 
depends sensitively on the local exchange couplings and 
thus provides a difficult challenge for the LDA. The dy¬ 
namical propagation of information and correlations in 
one-dimensional setups with cold atoms is a focal point 
of research at the moment and it is thus theoret¬ 

ically important to have accurate models for these sys¬ 
tems also in the strongly interacting regime where time- 
dependent exchange and non-equilibrium quantum mag¬ 
netism can be studied [5l|, . 

The paper is organized as follows. In Sec. El we out¬ 
line the model and its assumptions, and then discuss how 
to compute exchange coefficients exactly and within the 
LDA. Section|TTT]presents a comparison of exact and LDA 
results for double-well potentials. In Sec. IIVI we intro¬ 
duce the spin model picture of strongly interacting ID 
systems and we apply the local exchange coefficients in 
the context of quantum state transfer to investigate some 
consequences of the differences between different compu¬ 
tational schemes. Section |V] contains a summary, discus¬ 
sions, and outlook. 


II. FORMALISM 


We consider particles that are confined to move in one 
dimension (ID) and assume that along the direction of 
motion there is a trapping potential, V{x), which is the 
same for all the particles. The particles have short-range 
interactions that we model by a Dirac delta-function and 
in turn the Hamiltonian may be written in the following 


way 
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where pi and Xi are the momentum and coordinate oper¬ 
ators of particle i, m is the mass of a particle and V{xi) 
is the trapping potential for the zth particle. The inter¬ 
action strength is parametrized by gm- In experimental 
setups the ID confinement is achieved by applying a very 
tight transversal trap. It may then be shown that gio 
can be directly related to the scattering properties of the 
non-trapped atoms and is a function of the low-energy 
three-dimensional scatter ing length 0313 and the trans¬ 
verse trapping length a_L [1^ . What is very interesting is 
that one finds resonances where giD diverges due to the 
presence of the transverse confinement. This has been 
clearly demonstrated in recent experiments (^ . In par¬ 
ticular, the regime where \giD\ is very large is accessible 
experimentally (^ . For simplicity we will use the nota¬ 
tion g = giD from now on since this can give rise to no 
ambiguities in the present context. 

In the present paper we will consider two-component 
Fermi systems with N = Nf + atoms, where and 
Ni are number of atoms of components with spin projec¬ 
tion up and down, respectively. The zero-range two-body 
interaction in the Hamiltonian will act only between pairs 
with opposite spin projection. For pairs with the same 
spin projection the Pauli principle requires antisymmetry 
upon exchange of the two particles. As the wave func¬ 
tion must also be continuous, we may infer that when two 
particles with identical spin projections are close to each 
other the wave function has a continuous first deriva¬ 
tive. In turn, the Dirac delta-function two-body interac¬ 
tion has no effect. In experiment with atoms there are 
interactions between pairs of atoms with identical spin 
projection, but they are highly suppressed due to their 
short-range nature and thus can be safely neglected for 
our purposes. The general Hamiltonian above takes this 
implicitly into account as we ensure antisymmetry among 
like components in our iV-body wave functions. We note, 
however, that many of our results may be easily trans¬ 
ferred to two-component bosons with uniform interac¬ 
tions, i.e. a single g controlling the interactions between 
both identical and different components [s^. 

We now consider the strongly interacting regime where 
l/g —>■ 0 jH, m, [ 5 ^ . The most general eigenstate of the 
Hamiltonian in the limit 1/g —>■ 0 is 


afc6»(xp^(i),... ,xp^(Ar))^o(^i, ■ - ■:Xn), (2) 

k 

where the sum is over the N\ permutations, Pk, of 
the coordinates, Ofc are the coefficients that depend on 
the ordering of the particles, 0{xi,..., xn) = 1 , when 
xi < X 2 < ■ ■ ■ < xn and zero otherwise. 4'o is a fully 
antisymmetric Wparticle Slater determinant wave func¬ 
tion constructed from N single-particle wave functions 
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that are obtained by solving the corresponding single¬ 
particle Schrodinger equation with the potential V{x). 
In this paper we are interested in the lowest energy man¬ 
ifold of A^-body states which are obtained by taking the 
single-particle states with the N lowest energies. 

The wave function ^I/q describes the non-interacting 
iV-fermion system with the energy Eq. It is impor¬ 
tant to note that this iV-body energy is = 

N\/{N-flN^^l) times degenerate, and thus there is a mani¬ 
fold of M{N-f, N^) degenerate 7V-body states in the limit 
1/5 —?► 0 which is a quasi-degenerate manifold at large 
but finite g (which is where we will be working below 
as we consider quantum state transfer). This degeneracy 
arises from the fact that in the 1/g ^ 0 limit the parti¬ 
cles become essentially impenetrable. Yet, there are still 
M{N^, distinguishable ways that the particles may 
be ordered on a line which all have the same energy. In 
practice one may think of these various orde ring s of the 
spins along a line as a set of basis states (ssl . 1571. [58| . 

In the limit I/g —0 we can write the A^-particle en¬ 
ergy, E, to the linear order in I/g as [3^ 1^. 


E — Eq 
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where Aj = ~ important observation 

is that there is an Ofc-independent coefficient, a^, in this 
expression. It is a geometric factor that depends solely on 
the total number of particles N and on the potential V {x) 
and its single-particle eigenstates. Remarkably, it does 
not depend on what system one is considering as long as 
the interactions are strong, i.e. as long as we consider the 
limit 1/5 —>■ 0. If one is able to compute aj for a given 
V{x), then this may be used to study multi-component 
Fermi or Bose systems with all possible combinations of 
internal components among the particles. The geomet¬ 
ric factor Uj can in a certain sense be thought of as the 
local exchange coupling in the system. This interpreta¬ 
tion is very useful when mapping the system onto a spin 
model [H, [13, m, 113 ■ If one thinks of the N particles 
sitting on a line, then the index j on aj corresponds to 
a pair of particles and the exchange coupling for that 
pair is proportional to aj. Since strongly interacting ID 
systems are governed by exchange processes, the statics 
and dynamics of such systems is essentially dictated by 
the aj coefficients and we would thus like to compute 
these in as general circumstances as possible. Therefore, 
the coefficients aj will therefore be the main focus of our 
discussion. 

An exact expression for aj was first derived in (55| 


& f n*^i dx,6»(xi,..., xjv)S(xi - Xj) (9^'o)" 
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(Vk|d^) 


where the derivative in the integral denotes 
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where one first takes the derivative of the A^-body an¬ 
tisymmetric function dtp with respect to Xi and then 
subsequently sets xi = xn- The exact expression is an 
{N — I)-dimensional integral and the numerical calcula¬ 
tion of aj is by no means an easy task. It is therefore de¬ 
sirable to consider whether appropriate approximations 
can be made to access these quantities also for larger val¬ 
ues of N. An important observation was made in Ref. [13| 
where it was noticed that for N < 6, a, local density ap¬ 
proximation can be used in computing aj and this yields 
results that are off by only a few percent for the bench¬ 
mark case where V (x) is a simple harmonic oscillator 
potential. Using a highly accurate ansatz wave function 
for the Y-body problem, it later became possible to get 
a highly accurate approximation to aj for the harmonic 
trap for any N given as a ratio of quadratic polynomials 
in N [13. 

Here we are concerned with the question of how well 
the local density approximation (LDA) does for different 
potentials. We must therefore define and discuss how the 
LDA may be applied in the context of strongly interact¬ 
ing particles in ID. This discussion closely parallels that 
of Ref. [531 ■ The main inspiration for the LDA method 
in the present context comes from earlier work on the 
Hubbard model using the Bethe ansatz where Ogata and 
Shiba have shown that in the strongly interacting limit 
the spin and charge dynamics decouple [93 . The spin de¬ 
grees of freedom may correspondingly be described by a 
spin model of the Heisenberg type (we return to this later 
on) with an exchange coupli ng that is proportional to the 
third power of the density [ionl - [In3 . This result is de¬ 
rived for a homogeneous system with periodic boundary 
conditions which are the typical basic conditions needed 
to solve the Bethe ansatz equations [^. In order to 
transfer these results into the present context with non- 
homogeneous confinement of the ID system, Ref. [53| 
suggested to use the density from the LDA in the expres¬ 
sion for the exchange coupling in the spin model. This 
yields the expression 


(LDA) 
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Here UTpiXi) is the ID Thomas-Fermi density [13,163 
which is given by 


nprix) = -^\/2m(p - U(x)), (7) 

TTn 

where p is the chemical potential of the system. Ref. (53 | 
proposed to calculate the Thomas-Fermi density in the 
center-of-mass positions Xi of the *th and (* -I- l)th par¬ 
ticles. The position Xi can be written as [53| X^i = 
4 J dx(p(®)(x) -I- p^*+^^(x)), where 


pW (a;) 


J dxi... dxN0{xi ,..., xn)S{x — Xi)4'o 
(4^14') 
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dd/n 


dxi 


Taking this position makes sense in light of the inter¬ 
pretation of a as a local exchange coupling of a pair of 












4 


particles and thus by symmetry one should take the den¬ 
sity at their common center-of-mass which may be deter¬ 
mined within the iV-body system. One may argue that 
any discrepancy between the exact value of a and the 
LDA approximation can be alleviated by simply picking 
the right values of Xi. However, finding the right Xi 
appears as difficult as calculating the exact a directly. 
Note also that for large N finding is an equivalently 
difficult task and thus one needs some method of deter¬ 
mining densities accurately. Here we will ignore these 
problems and be concerned with how well the LDA per¬ 
forms in cases where we can compute both the exact and 
the LDA expression for a easily. We therefore work with 
N < 6 here. For a commonly used harmonic oscillator 
trapping potential V{x) = the Thomas-Fermi 

density can be written down in a simple form 

nTF{x) = —\ 2m{Nfiuj — —niLo'^x'^). ( 9 ) 

nn V 2 

In this case the approximation is rather accurate com¬ 
pared to the exact calculation and the relative error is 
not larger than 7 percent and decreases rapidly with in¬ 
creasing particle number N (^ . However, for potentials 
of more complicated shape the calculation of the center- 
of-mass positions becomes much more difficult. 

In the following we discuss the coefficients ai and 
^{LDA) ^ with * = 1,..., Ai — 1 for different potential pro¬ 
files. The numerical calculations of Eqs. (|3]) and ([5]) con¬ 
sist of two main parts. First, we construct the Slater 
determinant dto with the eigenstates, '0j(xi), of a par¬ 
ticle trapped in the potential V(xi). We assume that 
our system is confined in a box potential with length L. 
To handle this numerically, we use the normalized eigen¬ 
states 

i / \ . TTnXi 

(Pn(xi) = d— Sin (10) 

where n is an integer, as an expansion basis for the func¬ 
tions 

= ( 11 ) 

n 

To evaluate integrals in Eqs. o and Q we utilize a 
recursive many-dimensional integration procedure with 
the one-dimensional integrals being calculated using a 
standard trapezoidal integration routine. Even though 
more advanced integration methods of course are avail¬ 
able they are not readily available for the case of nested 
integrals (integration limits that depend on subsequent 
integrations). We find that the trapezoidal integration 
converges rapidly with the relative error for the values 
of the coefficients not exceeding 10“®. We hnd that 50 
basis vectors </>„ and 50 integration points are sufficient 
to achieve this accuracy in a reasonable time. In the 
following sections we dismiss the superscript [LDA) for 
simplicity and make sure that it is clear from the con¬ 
text of the discussion and in the figures whether we are 
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FIG. 1: The shapes of the potentials V\{x) from Eq. (I12II (top 
panel) and V 2 (x) from Eq. (1131) (bottom panel) for different 
values of parameters. 


discussing the exact expressions for a or the LDA ap¬ 
proximations. 


III. COMPARISON OF EXACT AND LDA 
COEFFICIENTS 

To perform a comparison between the LDA results and 
the exact solutions for the exchange coupling constants ai 
we use two specific forms of a double well potential which 
have flexibility to explore the similarities and differences 
between the LDA and the exact calculation. The first 
potential has the form 

Vi{x) = ^k{\x\ - b)^, (12) 

where k and b define the energy scale, and the height of 
the central barrier. Notice that this potential has a cusp 
(discontinuity of the first derivative) at a; = 0. This is 
perfectly well allowed in the formalism as this will not 
cause any troubling discontinuities in the single-particle 
wave functions or their first derivatives. The second po- 
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FIG. 2: Geometric exchange coupling coefficients ai for the 
double well potential in Eq. m with fc = 1 as function of 
the central barrier height b. The exact values based on the 
expression in Eq. @ are shown as solid fines, while the cor¬ 
responding approximate values based on the LDA formula in 
Eq. m are the dashed lines. The top panels shows N = 2,3 
and 4, the middle panel has N = 5, and the bottom panel has 
N = 6. Note that the parity symmetry of Eq. (O implies 
that ajv-i = cti. The dashed horizontal lines in the N = 4 
and N = 6 cases correspond to the value ai for N — 2 and 
N = 3, respectively. This is the excepted asymptotic values 
for large b as discussed in the text. 


tential is a symmetric trap of the form 


V 2 ix) = - Vo sin^ 


TT 

(^2a; A] 

. 2 


r2x A] 
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— usm 

TT 

vT + A 


(13) 

where the values of Vq and u control the shape of the 
potential. Figure [1] shows the shapes of the potentials 
we use in the article for different values of the control 
parameters. 

We note that we consider values in the interval x G 
[—L/2; L/2] where L is the length of a box which con¬ 
fines the whole system. This is required for our numer¬ 


ical procedure where we use L = dvr in our calculations. 
Energies are measured in units of e = /mL^, with a 

factor of 4 coming from the box extending to ±L/2 For 
the potential in Eq. (fT^ we set k = 1 (in units of L and 
e) while for Eq. m we have used Vq = 50 in units of 
e. In deep enough potentials the single-particle states do 
not feel the effect of this ’outside’ boundary. However, 
for shallow potentials some of the states might reside in 
the whole box and thus be influenced by the outer wall. 
We will see how it affects the exchange constants in the 
subsections below. 


A. Results 


We now consider the coefficients ai using Eqs. (|4|) and 
m for different values of the control parameters; the 
height of the central barrier b for the potential in Eq. m 
and the parameter u for the potential in Eq. ([0)1 . Figs.H] 
and |3] show both the LDA and exact values of the coef¬ 
ficients ai for different numbers of particles TV = 2 — 6. 
In this article we consider only spatially symmetric ID 
potentials and so a^-i = oti. 

In Fig. [2] we show the results for the potential in 
Eq. (foil comparing the exact results obtained from 
Eq. (HI) to the LDA results obtained with the formula 
in Eq. ([HI) for particle numbers TV = 2 — 6. Note that 
we do not need to plot all the ai coefficients as the par¬ 
ity invariance of the double well in Eq. m gives us the 
convenient symmetry relation ajsi-i = oti. The results 
demonstrate that the LDA does very well for small val¬ 
ues of the central barrier height b which should not be 
too surprising as the potential for small b resembles very 
much that of a harmonic oscillator which was previo usly 
shown to be quite accurately described within LDA . 

One notices that for the even particle numbers the mid¬ 
dle coefficient, i.e. ai, a2, and a^ for TV = 2, TV = 4, and 
TV = 6, respectively, has similar behavior, and that this 
behavior is not captured well by the LDA for b > 0.75. In 
fact, this middle coefficient, aAr/ 2 , calculated via Eq. ([6]) 
goes to zero much faster than the exact value. This q; 7 v /2 
pertains to the exchange coupling in the middle of the 
trap, i.e. to the exchange taking place right at the cen¬ 
tral barrier in the potential. It is therefore determined by 
the probability for particle tunneling through this central 
barrier, and we clearly see the LDA fail to capture this ef¬ 
fect for larger barriers, i.e. larger values of b, where LDA 
underestimates exchange. This is connected to the fact 
that LDA and in particular the Thomas-Fermi density 
in Eq. o is only well-defined in between the classical 
turning points, where /i — V{x) > 0. Hence for large 
barriers where fi — V(x) < 0, the Thomas-Fermi density 
goes to zero and thus in turn the coefficient ajv /2 goes to 
zero. For example, in the case of the potential in Eq. m 
= 0 for 6 > \/N (when fc = 1 as we have here). 

The termination values of 6 = \/2, 6 = 2, and b = ^/& 
can be clearly seen on the horizontal axis for ai (TV = 2) 

















6 



0 10 20 30 40 50 


u[e] 



u[e] 



u[e] 


FIG. 3: Geometric exchange coupling coefficients ai for the 
smooth double well potential in Eq. (US) as function of the 
parameter u that controls the central barrier. Here we have 
set Vb = 50. The exact values based on the expression in 
Eq. @ are shown as solid fines, while the corresponding ap¬ 
proximate values based on the LDA formula in Eq. are 
shown as dashed lines. The top panel shows N = 2,3 and 
4, the middle panel has N — 5, and the bottom panel has 
N = 6. Note that the parity symmetry of Eq. dn implies 
thsitj QiJV_ i — 


and a-i {N = 4) in the top panel of Fig. [21 and for 
{N = 6) in the bottom panel of Fig. |21 
When the height of the barrier is large (large b) each of 
the two wells can be approximated by the harmonic oscil¬ 
lator potential. We see that in this case the exact values 
of the coefficients becomes almost constant as function of 
the barrier height. For even number of particles we intu¬ 
itively expect that half of the particles will be located in 
each well. In the case with = 4 this would imply that 
ai should go to a constant for large b, reaching the value 
corresponding to the 6 = 0 case with N = 2 (two par¬ 
ticles in a single harmonic well). This is clearly seen in 
the top panel of Fig. |2] where the horizontal dashed line 


marks the latter value. Likewise, for N = 6, we expect a 
split into two three-body systems and thus that both ai 
and as should be approaching the a\ value for 6 = 0 and 
N = 3. Again this is seen very nicely in the bottom panel 
of Fig. 121 the first two coefficients for iV = 6 approach the 
dashed horizontal line corresponding to TV = 3 and 6 = 0 
(notice that there are different vertical scales on the three 
panels in Fig. HD. 

For the spatially symmetric potentials considered here 
with an odd number of particles, we cannot use the same 
logic of division of the particles into the two wells as 
the barrier grows large as for even particle numbers. The 
overlap across the barrier of the single-particle wave func¬ 
tions will not vanish for odd particle numbers. Hence, 
the coefficients ai which depend solely on these single¬ 
particle wave functions will also not vanish. We clearly 
see in the top and the middle panel of Fig. | 2 ] that for 
6 > 2 the LDA result does not do a good job in describ¬ 
ing the exchange couplings. Again, this is caused by the 
fast decrease of the Thomas-Fermi density on which the 
LDA result relies as the barrier increases. This decrease 
of density clearly takes place much faster than seen in 
the exact results. We thus see that odd-even effects can 
be considerable in the comparison of the LDA and the 
exact method. 

In order to further explore the case with odd particle 
numbers we have tried to apply a small tilt to the po¬ 
tential in order to explicitly break the parity invariance. 
This could for instance be done by an additional term 
that is linear in x in the potential (fT^ . When this is 
done for an odd number of particles the values of the 
exchange couplings approach the values of the couplings 
of a smaller system. For instance, for the system con¬ 
sisting of TV = 5 particles, two of the particles will be 
located in one of the wells, with the interacting coeffi¬ 
cient approaching the value of ai for two particles in a 
harmonic trap, and the other three will occupy the other 
well, with ai = a 2 approaching those of a three-particle 
system in a harmonic trap. At the same time we find 
the intuitive result, namely that the interaction between 
these two subsystems will approach zero as the barrier 
increases. We will not discuss the introduction of slight 
symmetry-breaking terms any further here. All in all, for 
the double well potential (TT^ the LDA approach provides 
a reliable approximation to the exact values of the coeffi¬ 
cients ai for small central barrier heights (small values of 
6 ) and then becomes poor for larger barrier heights. An 
exception is found for even particles numbers where the 
system splits into two equal size groups that can be de¬ 
scribed as particles in two separate harmonic wells. Here 
LDA does approach the exact results for the split system 
asymptotically. 

To further explore the robustness and/or failures of 
applying the LDA to our problem system, we now change 
the potential into a different kind of double well shape 
which has the form given in Eq. m- What one should 
notice here is that the potential in Eq. m is in a sense 
shallow, i.e. it does not increase to infinity around its 
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N 



2 

0.02487 

0.01328 (47%) 

3 

0.06963 

0.04485 (36%) 

4 

0.14921 

0.10630 (29%) 

5 

0.27355 

0.20798 (24%) 

6 

0.45260 

0.35911 (21%) 


( fT V A 

TABLE I: Comparison of the exact calculation, ' 

(first column), to that of the LDA, (second column), 

for a potential consisting of a flat box with infinite walls. No¬ 
tice that in this case the exchange coefficients are independent 
of i, i.e. Oi = a. The percentages in parenthesis in the sec¬ 
ond column give the deviations of the LDA from the exact 
results. Notice that the LDA always underestimates the ex¬ 
change coupling coefficients. 


edges. Therefore the box potential that we have as a 
hard-wall boundary around the system could be felt by 
the particles. This is in fact the case for larger particle 
numbers as we will now discuss. 

First we consider the cases with N = 2 and N = 3 as 
shown in the top panel of Fig. [31 Here we see that the 
qualitative behavior of the ai coefficients is similar to the 
potential of Eq. (fT31) (shown in the top panel of Fig. [3]) 
with the two-particle case shown a steady decrease with 
the barrier parameter, u, while the three-particle case 
first decreases and then has a slight increase again at 
large values of u. This is analogous to the behavior in 
the top panel of Fig. [3] We also notice that the LDA does 
a very good job for the two- and three-particle systems 
for most values of u in the top panel of Fig. [31 For the 
N = 3 case this should be compared and contrasted to 
the case in the top panel of Fig. [3] where the same N = 3 
case for the potential in Eq. m demonstrates that the 
LDA deviates significantly from the exact result for larger 
barrier heights (large b values in that case). 

For more particles {N > 4) the LDA results start to 
deviate significantly from the exact results as we see in 
all three panels in Fig. (3] We see that already for four 
particles the coefficient ai differs significantly in the LDA 
and exact calculations. This difference only increases for 
larger particle numbers. The reason for this discrepancy 
can be traced back to the form of the potential in Eq. m 
and its shallow nature as compared to that in Eq. m- 
For the parameters we have chosen, the potential does 
not become deep enough for five or six single-particle 
eigenstates to be confined completely inside of the po¬ 
tential profile. This can be rephrased by stating that not 
all the five or six lowest single-particle energy eigenvalues 
are in fact negative. The most energetic single-particle 
wave functions that we require to build the totally an¬ 
tisymmetric function, 'I'oj thus feeling the outside of 
the potential. In particular, they are influenced by the 
box potential that surrounds our system. While one may 
alleviate this problem by going to deeper potentials (in¬ 
creasing Vb in Eq. (fT31) ) we choose Vb = 50 in order for 


us to study quantum state transfer in the last part of our 
paper and make contact with recent results that utilize 
the same parameters [58| . 

To gain further insights into the influence of the box 
and the behavior of LDA in this respect we may con¬ 
sider just the box potential on its own. In Tab. U we 
present the results of applying both the exact formula 
(first column) and the LDA version (second column )to 
the box potential. We clearly see a large discrepancy for 
the particle numbers we have studied. The deviation of 
the LDA from the exact result is given in percentages in 
the parenthesis in the second column of Tab. HI By doing 
a crude fit to the percentages we find that the deviations 
scale approximately with the particle number as a power 
law We thus find a quite slow convergence and 

for smaller system sizes the deviations can be significant. 
We can see that as the single-particle wave functions are 
distributed over the whole confining infinite square well 
potential the local density approach becomes rather in¬ 
adequate. The large discrepancies in the values of the 
LDA and the exact results for the ai coefficients of the 
shallow potential in Eq. m for particle numbers N > 5 
can now be better understood. They are in large part the 
result of the fact that the outer box boundaries do be¬ 
come important for these particles numbers and that the 
LDA does a poor job of describing particles in a box with 
a flat bottom. We do see that the LDA approximation 
gives qualitatively similar results to the exact calculation 
even for TV > 4, but quantitatively the LDA may differ 
substantially. We have checked that as one increases the 
depth of the potential in Eq. m one does indeed see 
better agreement. However, the deviations we identified 
as common for both potential profiles within the LDA 
remain. 


IV. QUANTUM TRANSPORT PROPERTIES 


In order to understand some potential effects that 
could be implied by using the LDA instead of the exact 
solutions for the exchange couplings, we now consider 
a dynamical protocol that has recently been discussed 
in the context of strongly interacting systems in ID. As 
mentioned earlier, one may indeed map these systems 
into effective spin models (m,[ 57|. This can be done for 

l60l | and Bose systems 


both Fermi systems |55ll57 
To keep the discussion concise, we will mainly consider 
the example of the two-component Fermi system here. In 
that case, the spin mapping is into the famous Heisen¬ 
berg spin-i model whose Hamiltonian (up to a constant 
energy shift) is given by 


N-l 

■S^+\ (14) 

i=i 

where is a spin operator where cr^ = 

is the vector of the Pauli matrices. Here the 
nearest-neighbor interaction coefficients are related very 
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FIG. 4: Quantum state transfer fidelity, F, for an = 4 two- 
component Fermi system with Ji obtained from the potential 
in Eq. m for M = 0 (first row), u = 12.5 (second row), 
u = 25 (third row), and u = 50 (fourth row) where u is in 
units of e. The other potential parameter is Vo = 50 for all 
panels. The (blue) solid line corresponds to the Ji coefficients 
obtained from the exact calculation and the (brown) dashed 
line corresponds to the LDA approach. Time is measured in 
units of h/t. 


simply to the local exchange coefficients and the coupling 
constant as Ji = —-y- This is another justification for 
using the term ’local exchange coupling’ for the cti above, 
i.e. that under the spin mapping they appear as nearest- 
neighbor couplings in what is equivalent to a spin chain 
Hamiltonian. For Bose systems, one may have more in- 


FIG. 5: Same as Fig. |4] but for N = 5. 


volved spin models as the coefficients for x-, y-, and z- 
direction spin operators are not necessarily the same (see 
Ref. [53 for further details). Below, we will make one de¬ 
tour from the uniform Heisenberg spin model in Eq. m 
in order to consider the so-called XX model which is a 
special case of the Hamiltonian in Eq. (fni) where all the 
terms with z-component operators are eliminated. 

Spin-i chains have been proposed as media for quan¬ 
tum state transfer about a decade ago |l03l - ll05 |. The 
quantum state transfer in these chains essentially corre¬ 
sponds to flipping a single spin at one end of the chain 
and then dynamically evolving the state such that one 
may project onto the state at which the single spin flip 
has reached the other end for all subsequent times. The 
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probability of finding the flipped spin at the other end 
is known as the fidelity of the quantum state transfer 
and if it reaches unity at some subsequent time we say 
that the system supports perf ect q uantum state transfer 
|l04| . In the original proposal [l03l | a spin chain with con¬ 
stant Ji = J was considered and it was shown that the 
fidelity could not reach unity in this case. Later stud¬ 
ies demo nstrated th at if the coefficients are chosen as 
Ji = Joi/N{N — i) with Jo some overall constant, then 
perfect state transfer is possible (although only within 
the so-called XX model). One would then have an ideal 
communications channel for quantum information. How¬ 
ever, it turns out to be exceedingly difficult to produce a 
system that fulfills these requirements on Ji. In Ref. 
Volosniev et al. proposed strongly interacting ID atomic 
systems as a possible realization of perfect state transfer 
and found that the potential in Eq. 0 can give rise 
to perfect state transfer in the XX model for u = 12.5. 
We will return to this below. First we need to define 
the spin state space and the fidelity that we consider. 
For a spin chain with one impurity (a single spin that is 
flipped) and TV — 1 majorities (in other words TV.|- = 1 and 
iVj, = — 1) it is natural to use the basis of spin functions 

{itu •••;), im • • • i), li • • • m), i; • • • ut)}. in this 

basis we define the fidelity of the quantum state transfer 
which is a function of time and is given by 


F{t) 




(15) 


The fidelity can be straightforwardly interpreted as the 
Hamiltonian acting on the initial state on the right (with 
the single flipped spin on the left edge of the systems) 
and then projection on the final state on the left (with 
the single flipped spin on the right edge). In practice, one 
expands the initial state on the set of eigenstates of the 
Hamiltonian, then constructs the state at all later times, 
and then projects onto the final state which also has some 
expansion in terms of the eigenstates of the Hamiltonian. 
Notice that the overall time scale of the transfer depends 
on ai and on g. In the strict limit where 1/g = 0, there 
is no dynamics as the particles are completely impene¬ 
trable and all orderings are eigenstates (all of which are 
degenerate). However, in the more realistic case where g 
is large but finite, our effective spin models work to linear 
order in 1 /g and can thus be used to study dynamics in 
the strongly interacting regime. The timescale of trans¬ 
fer will then depend linearly on g. One may think of the 
state transfer process as a set of subsequent flips of pairs 
of spins along the chain. Each of these local exchanges 
depend on ai, i.e. if Ui is large this happens fast and 
vice versa. We thus see that large barriers will suppress 
the transfer as expected. However, the process can de¬ 
pend rather delicately on the actual values of the local 
exchanges as we will now demonstrate. For instance, in 
the exact results we may see slow suppression of the lo¬ 
cal exchange with barrier height as in Fig. [5] in regimes 
where LDA would give exactly zero (when the chemical 
potential is below the barrier). This is one source of error 
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FIG. 6: Same as Fig. |4] but for N = 6. 


in the LDA that could carry into a dynamical protocol 
like state transfer in a severe way. 

Here we plot the fidelity of Eq. as a function of 
time for different potentials of the form given in Eq. m 
for the cases with iV = 4 (Fig. [4]), 5 (Fig. [5]), and 6 
(Fig. [5]) particles, respectively. All the three cases we 
present show the characteristic oscillatory behavior of 
state transfer fidelities as function of time [^. How¬ 
ever, one does indeed notice that the oscillations become 
more prominent and more irregular as N increases (sys¬ 
tem size or chain length in the spin chain language), and 
also as the control parameter for the central barrier u 
is increased. Perfect state transfer is never achieved for 
these parameters in a Heisenberg model of the type in 
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Eq. (HID without an external ma gnetic fiel d. This is in 
accordance with previous results The calcu¬ 

lations demonstrate that even when there are only small 
differences in the exact and LDA results, the fidelity can 
show large variations particularly for longer time inter¬ 
vals. This is probably due to an accumulation of phase 
factors in the system over time that tends to drive the 
exact and LDA results apart. However, we do see some 
instances of very good agreement as for instance in the 
second row of Fig. 0] where u = 12.5 and Vq = 50 in 
Eq. (fOp . Even in this case, noticeable differences be¬ 
tween exact and LDA results are seen, but the overall 
agreement is very good. For the larger particle numbers 
in Figs. [5] and [6] we may even notice a clear tendency for 
the LDA results to predict large fidelity spikes that are 
either at different times as compared to similar spikes us¬ 
ing the exact result, or in some instances the LDA shows 
spikes where the exact results have none (see for instance 
the fourth row in Fig. [S]or the third row in Fig. [5D- We 
thus conclude that in most cases, the LDA results can 
provide large deviations from the exact results in a dy¬ 
namical process such as quantum state transfer. 

In closing this section, we want to consider how well the 
LDA results do in the case where perfect state transfer 
is achieved. This can be achieved with parameters corre¬ 
sponding to the second row in Fig. 0] in a setup where we 
can discard the interaction in the Hamiltonian in 

Eq. m and thus reduce the problem to the Heisenberg 
XX model. This can be realized by either using a tai¬ 
lored external magnetic field applied to the system or in 
specific models with strongly interacting two-component 
Bose systems In Fig. [7] we compare the fidelities 

computed using the exact and the LDA values of the 
exchange coupling coefficients. We see that the perfect 
transfer is indeed obtained for the values calculated via 
Eq. (0D. However, the fidelity of the spin transfer de¬ 
creases with time for the coefficients obtained via Eq. 

The relative differences in the ai exchange coefficients are 
no more than 2 percent when comparing the exact to the 
LDA method. However, we still see that it affects the 
quantum transfer properties on longer time scales signif¬ 
icantly. 


V. SUMMARY AND DISCUSSION 

We have considered strongly interacting two- 
component systems in one dimension held in place 
by an external confinement. Such systems have an 
effective Hamiltonian that can be completely specified 
by computing a set of local exchange coefficients which 
may be interpreted as nearest-neighbor spin exchange 
interactions when the system is mapped onto a spin 
model Hamiltonian of the Heisenberg type. Computing 
these exchange coefficients as accurately as possible is an 
important yet also computationally difficult task. In the 
present paper we have explored two different approaches. 
One is a ’brute force’ calculation of a multi-dimensional 



0 100 200 300 


time 

FIG. 7: Quantum state transfer fidelity, F, for an A = 4 two- 
component Fermi system with Ji obtained from the potential 
in Eq. for u = 12.5 and Vb = 50 (in units of e) for the case 
where the z-components in Eq. (ED can be discarded so that 
we are in the so-called Heisenberg XX spin model. In this 
case perfect state transfer can be achieved [s^l . The (blue) 
solid line corresponds to the Ji coefficients obtained from the 
exact calculation and the (brown) dashed line corresponds to 
the LDA approach. Time is measured in units of ft/e. 

integral which gives the exact result (but is prohibitive 
for larger particle numbers) and the other is an approach 
inspired by the local density approximation (LDA) that 
could in principle be used to reduce the computational 
complexity. While previous studies have shown that 
the local density approach can be accurate at the level 
of a few percent for the case of a harmonic oscillator 
potential, we explore more complicated geometries 
consisting of two instances of a double-well potential. 

Our findings demonstrate that while the LDA does 
rather well with potentials that resemble a harmonic os¬ 
cillator, the exchange couplings do not have the right 
qualitative and quantitative behavior in the LDA when 
there are significant barriers as is typical of a double¬ 
well potential. In particular, the LDA cannot capture 
the right quantum tunneling processes across such barri¬ 
ers and may thus leave out important effects. We have 
shown that for small systems this can lead to an under¬ 
estimation of exchange by the LDA for potentials which 
have a more interesting structure than the simple single¬ 
well harmonic oscillator. We expect this observation to 
have influence on experiments with small system sizes in 
ID optical lattices. There one typically also has an ex¬ 
ternal overall smooth potential (approximately harmonic 
in shape). Thus, the potential seen by the particles is 
the superposition of harmonic trap and optical lattice, 
i.e. a ’smiling’ lattice potential. This is a very struc¬ 
tured potential and in the strongly interacting regime 
one could be in dire straits with a simple LDA approach 
as compared to the exact exchange couplings. In order 
to test out the LDA in a concrete physical process, we 
considered quantum state transfer of single spin flips in 
systems with four, five, and six particles. Here we found 
that while the LDA performs reasonably for some double- 
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well realizations there is amplification of the deviations 
of the LDA compared to the exact results in the transfer 
fidelity that can be significant and lead to large errors in 
both maximum fidelity values and the specific times at 
which these are attained. 

The current study has concentrated on small particle 
numbers where very accurate results can be obtained for 
the local exchange coefficients using multi-dimensional 
integration so that a comparison between exact and LDA 
results is possible for arbitrary potentials. Incidentally, 
as we discussed in the introduction, the system sizes used 
here are also of great current experimental interest in cold 
atoms. One would, however, like to study how these re¬ 
sults scale to larger particle numbers. This most likely 
requires alternative approaches not only to the exact for¬ 


mula but also to the LDA formula. In the current im¬ 
plementation it depends on the density of the system 
which is not easy to compute for arbitrary potentials. 
One may thus pursue an agenda of finding an alternative 
LDA method that obtains the density by some other and 
computationally much simpler approach, and simultane¬ 
ously explore how to get a computational reduction of 
the exact formula. 
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